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ABSTRACT 

Context. According to previous investigations, the effect of diffusion in the stellar atmospheres and envelopes of hot white dwarfs and 
subdwarf B (sdB) stars strongly depends on the presence of weak winds with mass-loss rates M < 10~"M o /yr. 
Aims. As in most of these stars with luminosities L/L Q ^ 100, no wind signatures have been detected, the mass-loss rates are unknown. 
In the present paper mass-loss rates are predicted from the original theory of radiatively driven winds. 

Methods. The method of solution is modified so that the usual parametrization of the line force multipliers is not necessary. This is 
important especially for very thin winds. In addition we checked whether a one-component description is justified. As a consequence 
of various simplifications, the mass-loss rates are expected to be overestimated. 

Results. Results are presented for effective temperatures in the range 25000 K < T cff < 50000 K and for various metallicities between 
solar and Z/Z Q = 0.01. For (pre-) white dwarfs and sdB stars a stellar mass of M, = 0.5M o is assumed. For fixed values of r e fr, M,, 
and Z, the results predict decreasing mass-loss rates with increasing surface gravity and an increasing dependence of the mass-loss 
rates on the metallicity. For white dwarfs with logg > 7.0 no wind solution exists even if the metallicity would be solar. Winds with 
mass-loss rates around 10~ H to 10~ 10 M o /yr are predicted for the most luminous sdB stars with surface gravities of logg 5 5.5, if the 
metallicity is not significantly lower than solar. For lower values of M metals decouple from hydrogen and helium. 
Conclusions. If weak winds with M 5= 10~ 12 M Q /yr exist, the metals cannot be coupled to hydrogen and helium. This should lead 
to additional changes in the surface composition, which have not yet been taken into account in the diffusion calculations with and 
without mass-loss. A possible scenario is the existence of pure metallic winds with mass-loss rates of M ^ 10~ 16 M o /yr and with 
hydrostatic hydrogen and helium. 
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1. Introduction 

The abundance anomalies in several types of chemically pecu- 
liar stars are believed to be at least partially due to diffusion 
in the stellar atmosphere and envelope. The predictions of dif- 
fusion calculations strongly depend on the presence of mass- 
loss. In several papers (Unglaub & Bues [TP98l [2000l |200T| ), 
we investigated the combined effects of gravitational settling, 
radiative levitation, and weak winds on the chemical composi- 
tion of hot white dwarfs and subdwarf B (sdB) stars. According 
to the results, these effects can explain the decreasing abun- 
dances of helium and metals during the cooling process of white 
dwarfs on the upper cooling sequence (with effective tempera- 
tures r e ff ^ 50000 K). The mass-loss rates were estimated either 
from scaling laws or were a free parameter. The winds were as- 
sumed to be "chemically homogeneous". If M[ is the mass-loss 
rate of an element and £ its mass fraction in the photosphere, 
then this assumption states that M/ = fyM, where M is the total 
mass-loss rate. Such a wind prevents (if M ^ 10~ n M G /yr) or re- 
tards (if M < 10~ n M Q /yr) gravitational settling. The present pa- 
per investigates whether the assumed mass-loss rates are consis- 
tent with the predictions of the theory of radiation-driven winds 
and whether these winds may be chemically homogeneous. 

In thin winds the one-component description is question- 
able, because the momentum redistribution via Coulomb colli- 
sions between metals (which are preferably accelerated due to 
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the radiative flux) and hydrogen and helium may not be effective 
enough (e.g. Owocki & Puis 2002 ; Krticka et al. 2003). This may 
lead to selective winds in which the metals are expelled from the 
star, whereas hydrogen and helium are left behind (Babel 1995 
I19961 I. In contrast to chemically homogeneous winds, selective 
winds should directly change the surface composition. 

In the absence of any mass-loss and other disturbing pro- 
cesses (e.g. convective mixing), an equilibrium between grav- 
itational settling and radiative levitation should be expected. 
Because of saturation effects the radiative force on an element in 
the stellar atmosphere and envelope depends on its abundance. 
The radiative force decreases with increasing abundance, so an 
"equilibrium abundance" can be found, for which the radiative 
force balances the effect of gravitational settling. This theory is 
described in Chayer et al. (|1995alfl995bl > and references therein. 
Dreizler & Wolff ( fT999l > and Schuh et al. (120021 have incor- 
porated the diffusion theory into their NLTE stellar atmosphere 
code, so synthetic spectra have been calculated for an abundance 
stratification, which is given from the equilibrium condition be- 
tween gravitational settling and radiative levitation. However, 
for hot white dwarfs the quantitative agreement with observa- 
tional results is not satisfactory in many cases. The presence of 
weak winds may be a possible reason for these discrepancies. 

Winds have been detected in pre-white dwarfs with effec- 
tive temperatures in the range 30000 K < T eff ^ 150000 K 
and surface gravities 3.5 ^ logg ^ 6.0 (cgs units). These ob- 
jects are in a post-asymptotic giant branch stage of evolution, 
and they evolve with approximately constant luminosities close 
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to 1O 4 L to higher effective temperatures. As derived from ob- 
servations and theoretical calculations (see e.g. Kudritzki et al. 
[20 061 Pa uldrach et al. [20 041 Tin kler & Lamers [20021 Heral d et 
al. 120531 Koesterke et al. [19981 Koe sterke & Werner [19981 and 
the review of Kudritzki & Puis 2000 ), the mass-loss rates are be- 
tween 10~ 6 and lO~ 9 M /yr. The mass-loss rates of white dwarfs 
on the cooling sequence (if winds exist at all) are unknown. 

The sdB stars can be identified with models of the ex- 
treme horizontal branch (EHB) stars with masses of M, * 
O.5M (Heber [19861 Saffer et al. |1994| i, effective temperatures 
of 20000 K < T eff < 40000 K, and surface gravities of 5.0 < 
logg ^ 6.0. Helium is usually deficient with number ratios of 
10~ 4 < He/H ^ 0.1 (e.g. Edelmann et al. [20031 Lisker et al. 
2005 ). However, these helium abundances are by at least one or- 
der of magnitude more than would be expected from the equilib- 
rium condition between gravitational settling and radiative levi- 
tation (Michaud et al. |1989l l. A possible explanation for this dis- 
crepancy is the presence of weak winds that retard gravitational 
settling (Fontaine & Chayer [T9971 Unglaub & Bues [19981 . For 
mass-loss rates M m lO~ 13 M /yr, within the lifetimes of sdB 
stars near the EHB (w 10 8 yr), the helium abundances would 
gradually decrease from the solar value to He/H » 10" 4 . 

The helium deficiencies in sdB stars are accompanied by 
anomalies in the abundances of heavy elements. Some of the 
most recent results of spectral analyses are from Geier et al. 
d2008l O' Toole & He ber (|2006l Edelm ann et al. (120061 
Blanchette et al. (2006), and Chayer et al. d20061 >. The hotter 
sdBs with T e ff > 30000 K show strong deficiencies in many 
cases especially of light metals like Al, Mg, O, and Si by more 
than a factor of 100 in comparison to the solar abundances, 
whereas enrichments of elements heavier than iron by a factor 
near 100 are not unusual. The abundances of iron and nitrogen 
are usually close to the solar value. If the scenario with diffu- 
sion and chemically homogeneous winds were the correct ex- 
planation of these compositions, it should be possible to find a 
mass-loss rate for which the anomalies of all elements can be 
explained simultaneously. For the case with solar initial abun- 
dance, the calculations of Unglaub & Bues (2001) show that, in 
the case of chemically homogeneous winds, helium should al- 
ways be more deficient than at least the elements C, N, and O, 
which have been taken into account. No mass-loss rate exists 
that leads to deficiencies in C and O by more than a factor of 
100, whereas helium is only deficient by a factor of ten. This, 
however, is not an unusual composition in sdB stars (e.g. Heber 
etal.EoUOl). 

Pulsating and non-pulsating sdB stars with similar stellar pa- 
rameters coexist in the r e ff-log g diagram (Charpinet et al. 2006 ). 
The proposed pulsation mechanism is associated with a local en- 
hancement of iron (or other iron group elements) in a mass depth 
of about 10 _7 M», which is indeed expected from the equilibrium 
condition between gravitational settling and radiative levitation 
(Charpinet et al. 119971 Fontaine et al. 2003 ). However, accord- 
ing to the calculations of Chayer et al. (120041 ) and Fontaine et 
al. (|2006), a weak wind with only M — 6 * lO~ 15 M /yr would 
be sufficient to destroy the iron reservoir within a timescale of 
10 7 yr. For higher mass-loss rates, this should happen in even 
shorter time scales. Thus for M * lO _13 M /yr, which would be 
required to explain the helium abundances, this pulsation mech- 
anism would become questionable. 

Mass-loss has been detected in sdO stars that are more lumi- 
nous than sdBs (Hamann et al. |1981| Rauch[1993 ). For sdB stars, 
up to now there has been no observational proof for the existence 
of winds. From a quantitative analysis of the Ho- line profiles of 
40 sdB stars (Maxted et al. 2001), a comparison of synthetic 



NLTE Ha line profiles from static model atmospheres with the 
observations have revealed perfect matches for almost all stars. 
Only in the four most luminous sdBs have anomalous Ha lines 
with a small emission at the line centre been detected, which are 
possibly signatures of weak winds (Heber et al. 2003 ). For the 
case with T eS = 36000 K, logg = 5.5, and logL/L = 1.51, 
Vink (2004) found a similar behaviour of Ho- from a spectral 
analysis of Ha with his wind code, if the existence of a weak 
wind with M « lO~ n M /yr is assumed. 

For several effective temperatures in the range 25000 K < 
T e tf < 50000 K and for various metallicities, mass-loss rates will 
be predicted according to the theory of radiation-driven winds 
from Castor et al. (119751 henceforth CAK). As the usual scaling 
laws may be unreliable for thin winds, the method of solution 
has been changed slightly (see Sect. 2), so it is taken into account 
that the radiative force is limited for small wind optical depth pa- 
rameters and tends to some maximum value. To simplify the nu- 
merical method, later improvements in the CAK theory are ne- 
glected. As discussed in Sect. 6, it should be possible to derive at 
least an upper limit for the mass-loss rates. For each wind model, 
we checked whether a one-component description may be justi- 
fied (see Sect. 3). The results for r eff = 40000 and 50000 K and 
several metallicities, which have been obtained with line force 
multipliers according to Kudritzki (2002 ), are presented in Sect. 
4. The results for r eff = 25000, 30000, and 35000 K (Sect. 5) 
were obtained with force multipliers according to our own calcu- 
lations (see Sect. 2.4.2). The predicted mass-loss rates are com- 
pared with the results of Vink & Cassisi (120021 ). 

In Sects. 6.3 and 6.4, we discuss how the abundance anoma- 
lies in sdB stars could be explained. The arguments are also 
relevant for hot white dwarfs and chemically peculiar main se- 
quence stars, which are reviewed by Smith (1996 ). A subgroup 
are the HgMn stars. Because they are characterised by very low 
rotational velocities and weak or non-detectable magnetic fields, 
they are one of the best natural laboratories for studying the 
competing processes of gravitational settling and radiative dif- 
fusion (Vauclair & Vauclair 1982). The HgMn stars are charac- 
terised by enhancements of heavy metals (e.g. Hg, Mn, Pt, Sr, 
Ga), deficiencies of some light elements (e.g. He, Al, N), and 
isotopic anomalies of metals. Some of the most recent papers 
about these stars are from Zavala et al. (2007) and Adelman 
et al. (120061 ). Similar to sdB stars, weak winds with mass-loss 
rates of 10~ 14 to 10~ 12 M o /yr could lead to abundance anomalies 
that are different from the ones obtained from the equilibrium 
condition between gravitational settling and radiative levitation 
(Landstreet et al. |19981 ). However, for the typical stellar parame- 
ters 10500 K < T eff < 16000 K and logg * 4.0, these winds can 
hardly be chemically homogeneous, because the radiative force 
is too low (Babel |1996l ). 

In chemically peculiar main sequence stars, the abundance 
anomalies depend on stellar rotation (see e.g. Vauclair [2003), 
on the presence of magnetic fields (Turcotte 2003 ) and of con- 
vection zones. In hydrogen-rich hot white dwarfs and sdB stars, 
these effects do not seem to be the most important ones. In gen- 
eral both are slow rotators (Karl et al. 120051 Koester et al. |1998l 
Heber et al. [T9971 Heber & Edelmann |2004t . O'Toole (|2005l ) 
finds no correlation between magnetic field strengths in sdB stars 
and abundance anomalies. In some (pre-) white dwarfs, mag- 
netic fields have been detected (Jordan et al. 120051 120071 and 
references therein). However, diffusion is effective in all white 
dwarfs and not restricted to a subgroup of them. A thin su- 
perficial convection zone with a mass depth of about 10~ 12 M» 
may be present in sdB stars, but only for helium abundances 
He/H > 0.01 by number (Groth et al. [19851 
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2. Calculation of mass-loss rates 



the wind optical depth parameter can be written as 



CAK introduced the dimensionless optical depth parameter 



/ dv 



t = cr e v th p I — 



(1) 



where p is the density and dv/dr is the velocity gradient. Here, 
cr e is the electron scattering opacity with 



cr e = — cr T 
P 



(2) 



6.6524 * 10 25 cm 2 is the Thomson cross section for 



where <x x 

electrons and n e the electron number density. The ratio n e /p de- 
pends on the composition and on the degree of ionization. Both 
are assumed to be constant throughout the wind. In addition, the 
wind is assumed to be isothermal with T = T e ff. Thus the mean 
thermal velocity v& of hydrogen is 



Vth = 



2k R T, 



B 1 eft 



m H 



(3) 



where k-g, is the Boltzmann constant and wh the atomic mass of 
hydrogen. The line radiative acceleration can be written in terms 
of the line force multiplier M (f) and the luminosity L of the star 



1 L 

c Anr L 



M(t) 



(4) 



In the original CAK theory, the force multiplier is parametrized 
according to M(t) = kr a . With constant force multiplier pa- 
rameters k and a and if the radiative acceleration due to elec- 
tron scattering is neglected, it follows that M ~ L~(see e.g. Puis 
et al. 120001 ). Usually a is between about 0.5 and 0.7. However, 
the results of such scaling laws are questionable for the case of 
thin winds. This is because, according to the parametrized form, 
the force multiplier and thus the radiative force may become 
arbitrarily strong if the wind optical depth parameter is small 
enough. As a consequence, even for such compact stars as white 
dwarfs, for which in fact no wind solution exists, the scaling 
laws derived for more luminous stars predict small, but non-zero 
mass-loss rates. 

For t — > the log M (?) - log t is not a straight line as would 
follow from the parametrized form, but the relation flattens and 
the slope -a of this relation approaches zero (see Sects. 2.4 and 
5.1). For given excitation and ionization equilibrium, the force 
multiplier smoothly tends to a maximum value M max . To take 
this into account, Kudritzki (2002) developed a method with 
variable force multiplier parameters to calculate wind models 
for extremely metal-deficient hot stars. The methods used in the 
present paper as described in Sects. 2.1 and 2.2 allow the calcu- 
lation of the mass-loss rate without parametrization of the force 
multiplier. However, in contrast to Kudritzki's method, later im- 
provements in the CAK theory like the finite disk correction fac- 
tor (Pauldrach et al. 119861 Friend & Abbott [1986] ) and changes 
in ionization in the wind are neglected. 

For fixed excitation and ionization equilibrium, the force 
multiplier depends on the wind optical depth parameter alone 
(see Sect. 2.4.2). With the equation of continuity 



M = Anr pv 

and the definition 

i dv 
D = r 2 v— , 
dr 



(5) 



(6) 



t = cr e V th 



M 

AnD ' 



(7) 



so t and thus the force multiplier are a function of M and D. 
Then, as described below, the solution to the momentum equa- 
tion is straightforward. The inclusion of the improvements in the 
CAK theory would lead to additional dependencies, because the 
finite disk correction factor, as well as the ionization equilibrium, 
change with radius. This may complicate the method of solution. 
Moreover, in weak winds these effect are probably less important 
than the effect of the shadowing of the flux by the photospheric 
lines discussed by Babel ( 1996 ), which is not taken into account 
in the original CAK theory and thus in the present calculations. 
The consequences of the various simplifications and their rela- 
tive importance will be discussed in Sect. 6. 

The momentum equation for a stationary wind is 



dv 



dr 



\_dp_ 
p dr 



GM, 



(1-F e ) + : 



(8) 



where dp/dr is the gradient of the gas pressure, G the gravita- 
tional constant, and F e is the Eddington factor: 



Lcr e 



AncGM* 



(9) 



For sdB stars, T e is only about 0.01. The solution to the mo- 
mentum equation is especially simple, if the contribution of the 
gas pressure is completely (see Sect. 2.1) or partially (Sect. 2.2) 
neglected. 

2.1. Momentum equation without gas pressure 

If the contribution of the gas pressure to the momentum equation 
is neglected and if the force multiplier is given in parametrized 
form with constant parameters, then analytical formulae for the 
mass-loss rate and terminal velocity can be derived (see Sect. 
8.7.1 of Lamers & Cassinelli. 119991 ). The method is similar for 
non-parametrized force multipliers; but then the equations must 
be solved numerically. If in Eq. (8) dp/dr is neglected, Eq. (4) 
for g lac | inserted, and the resulting equation multiplied with r 2 , 
then the momentum equation with the definition of D can be 
written as 



H = D + GM t (1 - F e ) - —M(t) = . 

4nc 



(10) 



With the present assumptions, this equation only depends on M 
and D. All other quantities are fixed. For a given mass-loss rate, 
it has either two solutions, one or no solution for D. The mass- 
loss rate according to the CAK theory corresponds to the case 
with one solution. Then the critical point of CAK degenerates 
and every point in the flow becomes critical. This solution fulfils 
the critical point condition dH/d(dv/dr) = 0, which is equiva- 
lent to 



d JL = o. 

dD 

With Eq. (10), it follows that 
Lcr e dM (t) dt 



1 - 



Anc dt dD 



= 0. 



(ID 



(12) 



Equations (10) and (12) are solved numerically for M and for 
the critical value D c . From insertion of M and D c into Eq. (7) the 
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critical optical depth parameter t c follows, which for completely 
neglected gas pressure is constant throughout the flow. The so- 
lution v (r) is obtained by quadrature of Eq. (6) with D = D c , so 
that 



v(r) 



H-f(>4) 



(13) 



where vo is the velocity at r = For r — > oo and with vo <k v a 
the terminal velocity is 



2D C 



(14) 



This gives the usual ft - type velocity law with /3 — i. 
2.2. The shooting method 

To neglect the gas pressure completely is only justified in the 
supersonic region, where the wind essentially is driven by the 
radiative force. If one is interested in the solution v (r) of the mo- 
mentum equation in the region around the sonic point and in the 
subsonic region, the gas pressure must at least be partially taken 
into account. With the equation of continuity and the equation 
of state for a perfect gas, the gradient of the gas pressure can be 
written as 



dp / a dv 2a" 

— = -p\ + — 

dr \ v dr r 



Here a is the isothermal sound speed. With T = T e s, it is 



(15) 



(16) 



where ji is the mean particle mass. The first term in brackets 
on the righthand side of Eq. (15) corresponds to an acceleration 
caused by the gradient of the gas pressure in planar geometry. 
This is the contribution to the momentum equation, which is im- 
portant in the subsonic region near the photosphere of the star. 
The term 2a 2 /r is the acceleration due to the area expansion. 

In winds with temperatures T ~ T e ^, the gravitational accel- 
eration g(r) exceeds this curvature term by more than a factor of 
100 near the photosphere. As, however, g(r) decreases propor- 
tionaly to r~ 2 , whereas the curvature term decreases only with 
r _1 , a radius exists at which both accelerations are equal. This 
is denoted as the Parker point (Lamers & Cassinelli, 1999} and 
is located far outside the star at distances of some hundreds of 
stellar radii. In a wide region between the sonic point and the 
Parker point, the contribution of the gas pressure to the momen- 
tum equation is small. As we intend to solve the momentum 
equation by outward integration from the photosphere, at least 
the first part of dp/dr, which is important in the subsonic re- 
gion, must be taken into account. The curvature term, however, 
will be neglected. Then with Eqs. (8) and (15) it follows that 



dv ( a 
dr \ v 



2 (1 - r e) + gmd ■ 



(17) 



In the complete momentum equation the term 2a 2 /r should ap- 
pear on the righthand side. The omission of this term simplifies 
the solution topology and the numerical method (see Sect. 2.3). 
As the density in the photosphere is known from the hydrostatic 
equation, the velocity vo at the inner boundary (at r — R», which 
is assumed at a Rosseland mean optical depth r = 2/3) can be 



obtained for a trial value of M from the equation of continuity. 
With this initial value, Eq. (17) is integrated outwards. The high- 
est mass-loss rate, for which a solution can be extended from the 
photosphere to infinity, is the one according to the CAK theory. 
For all cases discussed in the present paper, the mass-loss rates 
have been calculated from both methods described in Sects. 2.1 
and 2.2. The results are very similar, because the differences do 
not exceed numerical uncertainties. 



2.3. Example: 05 main sequence star 

To illustrate the methods described in the previous sections, we 
now calculate a wind model for an 05 main sequence star using 
the same stellar parameters as CAK: 

T e fi = 49290 K , logg = 3.94, M t =6OM , 

and L = 9.66 * 1O 5 L , T e = 0.4, cr e = 0.325 cm 2 g-'. As in CAK 
for the force multiplier it is assumed that 



M( t )= 4(-^ f 

30 \Vl2 



-0.7 



(18) 



The factor 1/VT2 appears because CAK used the optical depth 
parameter for a carbon ion, so that their Vth is the thermal velocity 
of carbon. Although M (f) is given in parametrized form, we do 
not make use of the parametrization during the calculations. 

Solutions for M = 10~ 4 and 10~ 5 M o /yr fail at velocities 
somewhat higher than the sound speed. Thus these mass-loss 
rates are too high. The highest mass-loss rate for which the solu- 
tion can be extended outwards to infinity and the corresponding 
terminal velocity are 



M = 6.4 * 10" 



Mo 
yr 



1527^. 

s 



For comparison, CAK obtained M — 6.6 * lO~ 6 M /yr and v,*, = 
1515 km/s. If the finite disk correction factor were taken into ac- 
count and all other assumptions are unchanged, then Pauldrach 
et al. dl986l obtain M = 3.5 * 10" 6 M o /yr and Voo = 5123 km/s. 
After multiplication of Eq. (17) with r 2 , with Eq. (4) for g m d, it 
follows that with the present assumptions Eq. (17) depends on 
M, D, and v. In the limit v » a, it is identical with the momen- 
tum equation with completely neglected gas pressure. Thus for 
high velocities the solutions of Eqs. (10) and (17) approach each 
other, whereas in the inner regions, where the gas pressure is es- 
sential, the differences are great, as can be seen from comparing 
the solid and dashed lines in Fig. la. In addition to the critical 
solution, shallow solutions also exist. They are characterised by 
lower mass loss rates, e.g. M = 10" 6 or 10~ 7 M o /yr and terminal 
velocities below the CAK value. 

In Fig. lb the solutions D of Eq. (17) are shown as a function 
of the velocity. In the subsonic region one solution for D (and 
thus the velocity gradient) exists. In the supersonic region, how- 
ever, two solutions exist for each mass-loss rate: a shallow solu- 
tion with a low value of D and a steep one with a high value of D. 
For mass-loss rates higher than the critical value, the shallow and 
steep solution meet, so that for velocities somewhat higher than 
a, no solution of the momentum equation exists any more. The 
critical value of M is the highest mass-loss rate for which both 
the shallow and the steep solution can be extended outwards. 
These solutions approach each other and for v » a they con- 
verge to the critical solution of the momentum equation with 
completely neglected gas pressure (in this case D is constant). 
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Fig. 1. Results for T eff = 49290 K, logg = 3.94, M* = 60 M Q . a) 
Solutions v(r) (in ms _1 ) obtained with the shooting method for 
M = 6.4* 10" 6 M Q /yr (solid line; the cross indicates the location 
of the sonic point), M = 10~ 4 , 10~ 5 , 10~ 6 and lO~ 7 M /yr (dotted 
lines). The dashed line is the critical solution of the momentum 
equation with completely neglected gas pressure, b) Solutions 
D (in cm 3 s~ 2 ) of the momentum equation (17) as a function of 
the velocity for the mass-loss rates as in Fig. la. The horizontal 
dashed lines represents the critical value D c obtained from the 
solution of the momentum equation with completely neglected 
gas pressure. In both figures the curves are labelled with log M. 



For subcritical mass-loss rates, the shallow and the correspond- 
ing steep solution converge to different values of D for v » a. 

The solution topology of the complete momentum equa- 
tion with the curvature term taken into account is more com- 
plicated (Bjorkman 1995). It has an X-type topology in D and 
several critical points exist. Then the shallow solutions meet 
at the Parker point and cannot be continously (in D) extended 
outwards. For this reason they are ruled out by CAK. Then for 
only one mass-loss rate a smooth and continuous solution exists, 
which reaches from the photosphere to infinity. To find this solu- 
tion, it is necessary to switch from a shallow to a steep solution 
at the CAK critical point. As in the present case with neglected 
curvature term shallow solutions can also be extented from the 
photosphere to infinity, only an upper limit can be set for the 
mass-loss rate from the time-independent equations. From time- 
dependent calculations, however, Feldmeier & Shlosman (2002) 
suggest that shallow solutions evolve in time towards higher 
velocities and mass-loss rates. Thus their physical relevance is 
questionable. 

As the steep solutions do not exist in the subsonic region, the 
outward integration of the momentum equation from the pho- 



tosphere must start on a shallow solution. As for the critical 
mass-loss rate, the shallow and steep solution are identical in 
the limit v » a (which would not be the case if the curvature 
term were taken into account), the critical solution can be found 
without switching from the shallow to the steep solution. This 
simplifies the numerical method, especially in the case of weak 
winds, in which the wind optical parameter is in a range where 
the log M (?) - log t relation is curved. Due to the X-type solution 
topology of the complete momentum equation, the integration 
should start at the critical point in this case. This is the method 
used by CAK. To set the initial conditions at the CAK critical 
point, the knowledge of the slope of the log M (?) - log t rela- 
tion (which corresponds to the force multiplier parameter a) is 
required. Therefore this method leads to complications, if this 
relation is curved and thus a is not a constant. 

As can be seen from Fig. lb, a solution with a subcritical 
mass-loss rate and a terminal velocity higher than the CAK value 
would require a jump in D from a shallow onto a steep solution 
somewhere in the supersonic region (higher values of D lead to 
higher terminal velocities as can be seen from quadrature of Eq. 
(6)). To see if a discontinuity in D may be physically reasonable, 
the time-dependent equations should be considered. As, how- 
ever, the CAK theory is based on the use of the Sobolev approx- 
imation for calculating the radiative force, the mathematical na- 
ture of the momentum equation would change if the Sobolev ap- 
proximation were dropped (Lucy 2007a; 2007b). As the Sobolev 
approximation is questionable at least in the subsonic region, it 
cannot be excluded that smooth and continuous solutions with 
subcritical mass-loss rates and high terminal velocities indeed 
exist in this case. 

2.4. The force multiplier 

The results for r eff = 40000 K and 50000 K, which will be pre- 
sented in Sect. 4, were obtained with force multipliers according 
to Kudritzki (2002). To avoid the methods of solution described 
in Sects. 2.1 and 2.2 being complicated by changing ionization 
and excitation equilibrium, the dependence of Kudritzki's force 
multipliers on the degree of ionization is only approximately 
taken into account. How this is done, is described in Sect. 2.4.1. 

The results for r eff = 25000, 30000, and 35000 K, which will 
be presented in Sect. 5, have been obtained with force multipli- 
ers from own calculations. These calculations and the various 
assumptions are described in Sect. 2.4.2. 

2.4.1. Force multipliers according to Kudritzki |2002) 

For r eff = 40000, 50000, and 60000 K, Kudritzki d2002b calcu- 
lated force multipliers with a line list of 2.5 * 10 6 lines of 150 
ionic species and with approximate non-LTE occupation num- 
bers as a function of the wind optical depth parameter and of 
n e /W, where n e is the electron density and W the geometrical 
dilution factor (W « 0.5 near the photosphere). In the present 
calculations n e /W is assumed to be fixed throughout the wind. 
The dependence on the degree of ionization is taken into account 
with the following iteration procedure. A first wind model is cal- 
culated with an input value (e.g. n e /W = 10 n cirT 3 ). From this 
model we obtain the mean value of n e /W at the sonic point and 
at a radius r = 2R t . The iteration is complete, if this output value 
is consistent with the input value within an accuracy of about 10 
%. Otherwise the iteration is repeated with an improved input 
value to ensure that the value of n e /W used in the calculation 
of the force multipliers has a similar order of magnitude to the 
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Fig. 2. Mass-loss rates (in M Q /yr) as a function of surface grav- 
ity for r eff = 40000 K, 50000 K, and Z/ Z = 1 calculated with 
force multipliers according to Kudritzki (120021 1. The solid lines 
represent the results obtained with the iteration procedure for 
n e /W as described in the text. The short-dashed and long-dashed 
lines represent the results for n e /W = 10 14 and 10 10 cirT 3 , re- 
spectively. 



value of the corresponding wind model near the sonic point and 
the CAK critical point (which according to the original theory 
with the complete momentum equation is near r = 1.5.R*). 

In Fig. 2 for T eS = 40000 and 50000 K and for Z/Z Q = 1, the 
mass-loss rates calculated with this iteration procedure are com- 
pared to the results obtained with the assumption that it is always 
n e /W - 10 I4 cirT 3 or 10 10 crrr 3 , respectively, regardless how 
high the wind density is. It can be seen that, for T e s = 40000 K 
and logg S 5.0, the results are approximately the same in each 
case with deviations less than a factor of two. This is because the 
dependence of the force multiplier on n e /W is weak for small 
wind optical depth parameters, which are expected in thin winds 
(see Fig. 3). The situation is different for T e g = 50000 K. Here 
the mass-loss rates derived for the two values of n e /W may dif- 
fer by almost two orders of magnitude. The values obtained with 
the iteration procedure are in between these results (the finally 
adopted values of n e /W can be read off from Fig. 5 in Sect. 4). 

For T e ff = 60000 K, the dependence of the mass-loss rates on 
the degree of ionization is even stronger; e.g., for \ogg = 6.0 it 
follows that log M = -1 1.4 with n e /W = 10 10 cm" 3 and logM = 
-7.6 with nJW = 10 14 cnT 3 . Thus the 

predicted value of M depends almost linearily on n e /W. 
Then the iteration procedure for n e /W does not converge. For 
this reason the case r e ff = 60000 K is not considered in the 
present paper. 



Fig. 3. Line force multipliers for T eff = 40000 K, Z/Z = 1 
as a function of the wind optical depth parameter according 
to own calculations (solid line) and according to Eq. (25) of 
Kudritzki d2002b for nJW = 10 14 (dashed line) and 10 10 cm 3 
(dotted line). Kudritzki's values have been not allowed to exceed 
log M max =3.37. 



The force multipliers according to Kudritzki's Eq. (25) have 
been not allowed to exceed a value of M max = 2000 Z/Z Q + 
Mn,He, according to his Eq. (12). The contribution Mu,n e of hy- 
drogen and helium to M max was estimated from own calculations 
as described below. The adopted values are Me.He = 263 for 
r eff = 40000K and M H , He = 93 for T eff = 50000K. At least for 
T e ff = 40000K, Z/Z = 1 this restriction is in good agreement 
with the value M max = 2340 obtained from own calculations (see 
Fig. 3). It only affects the results for extremely weak winds with 
wind optical depth parameters log t % -7.0. In these cases de- 
coupling of the ions from hydrogen and helium is expected, so 
the corresponding wind models are of questionable physical rel- 
evance. Thus this restriction does not change the conclusion of 
the present paper. 

2.4.2. Force multipliers from own calculations 

The results presented in Sect. 5 are obtained with force multipli- 
ers from own calculations. The elements H, He, C, N, and O are 
taken into account with a line list that is essentially similar to the 
one used in the diffusion calculations of Vauclair et al. ( 1 19791 1 
and later in our diffusion calculations (Unglaub & Bues, 1996). 
It consists of about 150 preferably strong lines with atomic data 
from Wiese et al. (1966). Similar to the diffusion calculations, 
multiplets, or doublets are lumped together into one line (note, 
however, that this may underestimate the radiative force, if not 
all lines are optically thin). All ionization stages are taken into 
account, with the exception of the neutral and singly ionized 
stage of the CNO elements. The occupation numbers are cal- 
culated with the assumption of LTE, for a temperature T = 
and an electron density « e = 2.0 * 10 14 cirT 3 , which is typical 
of regions near the wind base. As the typical electron densities 
in the outer parts of weak winds are considerably lower and be- 
cause the assumption of LTE is unreliable in the wind region, 
an agreement of the derived mass-loss rates with the results of 
more sophisticated calculations can be expected only if the de- 
pendence of the radiative acceleration on the ionization and ex- 
citation equilibrium is weak enough. 
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As in Abbott ( 119821 1 the force multiplier can be written as 
A/(/)= > -Am- — ' ' . (19) 



the sum is over all lines. If vo is the line centre frequency, the 
Doppler width is defined here according to 



Av D = — vth ■ 

c 



(20) 



Note that v t h is the thermal velocity of hydrogen defined in Eq. 
(3). With expressions (1) and (22) for t and rji, it can be seen that 
in Eq. (19) v± cancels out. 

The ratio F v /F indicates the monochromatic to total flux at 
the frequency of the line. This flux is calculated from similar 
models as used in the diffusion calculations e.g. in Unglaub & 
Bues (2001). The temperature structure is obtained with the as- 
sumption of LTE and with the diffusion approximation for the 
stellar flux: 



(21) 



with = iTgL at the outer boundary. Here, f is the Rosseland 
mean optical depth. The Rosseland mean opacity is calculated 
from the monochromatic continuum opacities for 264 wave- 
lengths in the range 0.1 A < A < 10000A. Similar to the calcu- 
lation of the force multipliers, the elements H, He, C, N, and O 
are taken into account, while F v is the emergent monocromatic 
flux of these models at the frequency of the line. It is important 
to note that only the continuum opacity is taken into account and 
not the opacity due to lines. Thus it is assumed that each line in 
the wind "sees" the continuum flux, which is independent of the 
velocity and thus the Doppler shift. Although this assumption is 
usual in the CAK theory, it is crucial for the case of thin winds 
(see Sect. 6.2). 

If stimulated emissions are neglected, then r\i is 



lie 



-f 



m e c pcr e Av D 



(22) 



where e is the electron charge, m e the electron mass, / the oscil- 
lator strength of the line, and n, the number density of particles 
in the lower level. With Eq. (2) for <x e , it can be seen that, for 
fixed degree of ionization and occupation numbers, the quantity 
77/ for each line is constant throughout an isothermal wind. As in 
addition the flux F v is fixed for each line, with these assumptions 
the force multiplier depends only on t. 

For t — > the force multiplier converges to its maximum 
value 



— Av D 77/ . 

lines 



(23) 



With expression (22) for rji it can be seen that, in this optical thin 
limit, the strong lines of the most abundant ions preferably con- 
tribute to the radiative acceleration. According to Abbott ( 1982) 
for small wind optical depth parameters log t S> -6 about 90 % of 
the total force multiplier originates from only about 100 strong 
lines. Thus for the calculation of extremely weak winds the small 
number of lines taken into account in the own force multiplier 
calculations may be sufficient. For higher wind densities, how- 
ever, the radiative acceleration and thus the mass-loss rate are 
underestimated, because in the optically thick limit the contribu- 
tion of a line is independent of its oscillator strength (e.g. Puis et 



al.|2000). Then the number of lines taken into account is espe- 
cially important. 

In Fig. 3 for r e ff = 40000 K and solar metallicity, the force 
multipliers from own calculations are compared with the ones 
according to Eq. (25) of Kudritzki d2002b for nJW = 10 14 and 
10 °cm _3 , respectively. For the higher values of f, the line list 
and the number of elements taken into account in the own calcu- 
lations is not sufficient. Thus the corresponding force multipliers 
are smaller than Kudritzki's ones. For low values of t the results 
are in good agreement. This shows that the force multipliers ob- 
tained with the present assumptions may agree with the results 
from more sophisticated calculations if the wind optical depth 
parameter is low enough. However, as the assumption of LTE in 
the present calculations may be unreliable, this agreement also 
requires the dependence on the ionization and excitation equi- 
librium to be small. This is obviously the case in the present 
example, because the differences between Kudritzki's results for 
n e /W = 10 14 and 10 10 crrT 3 decrease to lower values of t. 

3. The decoupling of elements 

According to our calculations, the mean radiative acceleration of 
the metals is by at least two orders of magnitudes higher than of 
hydrogen and helium. This may justify lumping hydrogen and 
helium together into one "element" (henceforth element 1) and 
the metals into "element" 2, although the velocities of the indi- 
vidual metals are different (Krticka 2006]). 

A chemically homogeneous, radiatively driven wind may 
only exist if the coupling due to Coulomb collisions between 
metals and hydrogen and helium is efficient. If a wind model is 
calculated from a one-component description it is implicitely as- 
sumed that this is the case. Whether this assumption is justified 
or not may be checked by considering the momentum equations 
of the individual constituents. 



3.1. Equations for an isothermal three-component model 

For each of the constituents, a momentum equation must be 
valid: 

dv\ dpi GM„ m . „ 

Pi y l 37 + -J7 + Pl - n ^ eE ~ PlSrad = A 2l < 24 ) 



dr dr r 1 - 
dv2 dpi GM„ 



P2V2 —j- + -rr- + P2 — 2 n 2 Z 2 eE - p 2 g ( ^ d = AQ 2 . (25) 

Here, pi, p 2 , Hi, n 2 , and Vi, v 2 are the mass densities, particle 
densities, and mean velocities of "elements" 1 (H and He) and 
2 (metals), respectively, and g ^ and g ® are the mean radiative 
accelerations and Z\e, Z 2 e the mean charges. And AQi and AQ 2 
are the momentum per unit volume and unit time, which is ex- 
changed between the two constituents via collisions. Because of 
momentum conservation, it is Agi + &Q 2 = 0. The Coulomb in- 
teraction of the elements with the electrons is taken into account 
via the polarization electric field E, which is obtained from the 
momentum equation for the electrons: 



dv e dp e 



dr 



dr r L 



+ n e eE - p e 



'rad 







(26) 



where g^y. is the radiative acceleration on the electrons due to 
Thomson scattering. 

The velocities of the various constituents are related to the 
centre of mass velocity v according to 

P1V1 +p 2 V'2 +p e v e 
v = . (27) 
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In addition we demand that there is no net electric current so that 
the plasma remains neutral on a macroscopic scale: 



Z\en\V\ + Zien2Vi — en e v e — . 



(28) 



This system of equations, together with the equations of conti- 
nuity for each constituent, and the energy equations (if the wind 
is not assumed to be isothermal) have to be solved to obtain 
a multicomponent wind model. This has been the approach of 
Krticka & Kubat (2000 2001a). In the present paper for each 
wind model obtained from the one-component description (from 
which M and v (r) is known), we check whether the momentum 
transferred from the metals to hydrogen and helium can be suf- 
ficient (see Sect. 3.3). Possible heating processes (e.g. frictional 
heating) are neglected, we assume always T = r e jf . 

3.2. The collisional acceleration on H and He 

A mean collisional acceleration on hydrogen and helium due to 
collisions with metals may be defined according to 



'coll 



1 

Pi 



-AQi 



From the equations derived by Burgers ( 1969) it follows that 

g^ ]1 =P2h 2 G(x) (30) 
with 



kn = 



1 \nZ\Z\e A 



m\m2 k&T 



In A. 



(31) 



With the assumed number ratio He/H = 0.1, the mean mass of H 
and He is m\ — 1.27m p . The metals are represented by the CNO 
elements with abundances proportional to solar ones (the solar 
abundances are from Grevesse & Sauval 1998). For the mean 
mass of the metals, m.2 = 14.6ra p is assumed. The Coulomb log- 
arithm In A is defined according to Burgers (119691) : 



1 (3k B TR D 
In A = — + In — 

2 \ Z { Z 2 e 2 



(32) 



Thus, in an isothermal wind with constant ionization, kn is al- 
most constant apart from a weak dependence on density via the 
Debye radius Rp, which appears in the Coulomb logarithm: 



fir 



k*T 



4ne 2 (riiZ 1 + n 2 Z 2 + « e l 



(33) 



As the particle densities decrease in an outward direction, the 
Debye radius and thus In A increase. The Chandraskhar function 



G W = -L(erf(x)-xA e xp(-x 2 )) 



depends on the difference in the mean velocities of the con- 
stituents. For v'2 > v'i it is 



x - 



with 



v 2 - Vi 



a 



2k B T 



(35) 



(36) 



where a has the dimension of a velocity and should not 
be confused with the force multiplier parameter, and ntu = 



m\m2l (mi + mi) is the reduced mass. As m\j — 1.17ra p in the 
present case, the quantity x approximately is the velocity differ- 
ence between both constituents in units of the thermal velocity 
v t h of hydrogen. The error function erf (x) has been evaluated 
with the routines in Press et al. ( 119921 ). For x <sc 1 the func- 
tion G (x) increases linearly with x. For x = 0.968, however, it 
reaches its maximum value 

G max = 0.214 (37) 
and decreases to higher values of x. 

3.3. Criterion for the existence of a coupled wind 

In a chemically homogeneous wind with a total mass-loss rate 
M, the mass-loss rates of the various constituents are by def- 
inition (see Sect. 1) Mi = ([M for hydrogen and helium and 
M 2 = ( 2 M for the metals. Then the equations of continuity for 
both constituents can be written as 



(29) ( X M = Ajir 1 p l v l 



£ 2 M = 4nr z p 2 v 2 ■ 



(38) 



(39) 



The mass fractions £\ and £2 of hydrogen and helium and metals 
specify the metallicity of the star. For solar composition we use 
(2 = 0.0133. As M is the total mass-loss rate, these continuity 
equations state that the relative mass flows of the constituents 
correspond to their relative mass fractions in the stellar atmo- 
sphere. Only if this is true can the surface composition of the 
star be expected to remain unchanged. 

In the present paper, metals are considered as trace elements 
so that p2 <K p\. As it is p e <sz pi, it follows that v\ w v from 
Eq. (27). The velocity v>2 of the mean metal must be higher than 
vi for that momentum is given from the metals to hydrogen and 
helium. Thus, with V2 = v, Eq. (39) yields an upper limit for the 
mass density p2- With this value and with G(x) = G max , from 
Eq. (30) an upper limit for the mean collisional acceleration on 
hydrogen and helium can be derived: 



(max) _ &M 



(40) 



For each wind model obtained from a one-component descrip- 
tion and with v\ » v and p\ ~ p the lefthand side of the momen- 
tum equation for hydrogen and helium (Eq. 24) is known at each 
grid point. Then the constituents may be coupled only if 



(max) 



(34) Scon 



dv 1 dpi 
dr p dr 



GM t 



1 

m\ 



-ZieE-g 



(1) 

rati 



(41) 



If this condition is not fulfiled eyerywhere in the wind, then the 
wind solution is incompatible with the assumption that the wind 
is chemically homogeneous. In this case clearly a multicompo- 
nent description would be required. The mean radiative acceler- 
ation on hydrogen and helium is obtained from own calcu- 
lations with assumptions as described in Sect. 2.4.2. The results 
have shown that its effect is negligibly small. In Sects. 4 and 5, 
mass-loss rates are predicted as a function of surface gravity. If 
#rad were neglected in criterion (41), the maximal surface grav- 
ity up to which coupled winds can exist would be lower by about 
0.05 dex only in logg. 
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3.4. Comparison with criteria from other authors 

In condition (41) the contribution of the gas pressure to the mo- 
mentum equation is taken into account. In the supersonic region, 
where this contribution is small, criterion (41) reduces to the one 
derived e.g. by Owocki & Puis (2002). If in Eq. (26) all terms 
proportional to p e are neglected, it follows that 



E « - 



1 d£e 
en, dr 



(42) 



With this approximation the radiative acceleration g due to 
the light scattering on free electrons has been neglected. This 
is justified in stars with F e « 1 (see Sect. 2). If the metals are 
trace elements, it is rc e as Z\ti\. With Vi ~ v, from Eq. (28), 
it follows that v e ~ v. Thus hydrogen and helium, as well as 
the free electrons, move approximately with the centre of mass 
velocity v. If approximation (42) is inserted into condition (41), 
then with p ~ p\ + p s the latter can be written as 



(max) . dv 

? „ > V 

'coll - dr 



GM, 1 dp 
r 2 p dr 



'rad ' 



(43) 



If the contributions of the gas pressure and of the radiative force 
on hydrogen and helium are neglected, this criterion is equiva- 
lent to the one in Owocki & Puis (2002). This is a good approxi- 
mation in the supersonic region, whereas the criterion used in the 
present paper may also be applied in the subsonic region, where 
the contribution of the gas pressure is essential. If the accelera- 
tion term vdv/dr is also neglected, it follows that 



(max) GM t 

<?coii - ji 



(44) 



In the supersonic region, this is a necessary but not sufficient 
condition for the existence of a coupled wind. It only requires 
that hydrogen and helium do not decelerate. In wind models ob- 
tained from the improved version of the CAK theory, the accel- 
eration term should be more important than the gravitational one 
(Gayley 2000). In the present wind models for high gravity stars 
obtained from the original version of the CAK theory, both terms 
have a similar order of magnitude in the supersonic region. 

For a given mass-loss rate and metallicity, with Eq. (40) for 
8 iT ' a max i mum velocity can be obtained up to which criterion 
(44) may be fulfiled: 



Vmax — ~. , , . . ^-12G max ■ 
47TGM, 



(45) 



This v max is an upper limit of the velocity, up to which the 
constituents may be coupled. In accelerating winds decoupling 
should in fact occur at lower velocities. Equation (45) corre- 
sponds to Eq. (22) of Owocki & Puis (2002) for the case w = 
in their notation. 

This maximum velocity may be compared with the surface 
escape velocity: 



2GM, (1 - T e ) 



R 



For T e «c 1, with R t = ^GM t g~ l , it follows that 
v esc = V2(GM,g)* . 
With the definition 

«2 

y = — 

"i 



(46) 



(47) 



(48) 
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Fig. 4. Solutions v (r) (inms 4 ) for T eS = 40000 K, logg = 5.8, 
M, = O.5M , and (from top to the bottom) Z/Z = 1,1/3, 0.2, 
and 0.1. Dashed lines indicate decoupling of metals, the dotted 
line represents the surface escape velocity. 

and (2 ~ m2in~ l y, with Eq. (31) for k\2, the ratio v max /v esc can 
be written as 

Vmax „> Z\z\e A \aKG u 
= My- 

Vesc 



m 2 k B T 



V2 



(49) 



For the wind models considered here, the Coulomb logarithm 
InA varies between about 6.0 at the wind base and 18.0 in the 
outer regions of the thinnest winds. With the values assumed by 



Owocki & Puis (2002), In A = 20, Zi = 1 , m x 
mass M» = O.5M as considered here, it is 



— = 1.4* 10 20 MyZiT- L g 



72 T -i„-i 



m„, for a stellar 



(50) 



where the mass-loss rate M is in M /yr, and the surface grav- 
ity g in cms 2 . The temperature T (in K) in the present paper 
is assumed to be equal to the effective temperature of the star. 
Frictional heating could reduce v max and thus favour decoupling, 
if the effect of a higher temperature is not compensated by a 
higher mean charge of the metals. 



3.5. Example: r eff = 40000 K, logg = 5.8 

For r eff = 40000 K, logg = 5.8, M, = 0.5M o , and several metal 
abundances, the solutions v (r) are shown in Fig. 4. It is indicated 
in which parts of the solution condition (41) is fulfiled, and the 
maximum velocities va up to which this is the case are given in 
Table 1 . In these calculations the metals are represented by the 
CNO elements with force multipliers from own calculations as 
described in Sect. 2.4.2. 

The only case where the constituents may be coupled 
throughout the wind is for Z/Z = 1, for which a mass-loss rate 
logM — —11.3 is predicted. (This value agrees with the one ob- 
tained in Sect. 4 with the force multipliers from Kudritzki 2002). 
This mass-loss rate is lower by a factor of ~ 1000 than in the 
main sequence star rSco, for which Springmann & Pauldrach 
(fT992l used the values M = 5 * lO~ 9 M /yr, M, = 19.6M , 
R* = 5.5R Q and T e ff = 33000 K in their analysis of multicompo- 
nent effects, which corresponds to logg = 4.3. As in the present 
example the stellar radius is only R„ = Q.15R Q , the mass flux 
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Fig. 5. Left: Predicted mass-loss rates (in M /yr), terminal velocities, and the adopted values of n e /W (in cm ) as a function of 
the surface gravity for T eff = 40000 K, M. = 0.5M Q , and Z/Z = 1 and 0.2, respectively. Right: The same for 7/ eff = 50000 K and 
Z/Z = 1, 0.2, and 0.01. The results for which decoupling of the metals is expected are represented by dashed lines. 



Z/Zo 


logM 


v„/v esc 


Vd/Vesc (41) 


Vraax/V csc (50) 


1.0 


-11.3 


1.10 




5.64 


1/3 


-11.9 


0.63 


0.186 


0.47 


1/5 


-12.9 


0.35 


0.017 


0.028 


1/10 


-14.7 


0.26 


0.0004 


0.0002 



Table 1. Predicted mass-loss rates (in M /yr), terminal veloc- 
ities, and velocities at which decoupling is expected accord- 
ing to criteria (41) and (50), respectively, for T eB = 40000 K, 
logg = 5.8, and various metal abundances. The surface escape 
velocity is v esc = 1136km/s. 



Ml (4ttRI) near the stellar surface has a similar order of mag- 
nitude to the one in tSco. Thus the wind densities are similar. 
Springmann & Pauldrach ( [19921 and Krticka & Kubat d2001bl > 
find that in rSco the multicomponent effects are important due 
to their contribution to the energy equation. In the present exam- 
ple, the relative velocities between metals and the passive plasma 
are up to 7 km/s. This may lead to frictional heating. Thus, as 



in rSco, the present example is near the border of the runaway 
condition. Multicomponent effects may have some importance; 
however, the coupling should still be effective enough the pas- 
sive plasma to be expelled from the star. 

If the metal abundance is reduced by a factor of three, 
log M = -1 1.9 is predicted by the one-component wind model, 
which is lower by a factor of four as in the case with Z/Z = 1. 
Because of the reduced mass-loss rate and the additional reduc- 
tion of the metal abundance, the density of metals in the wind 
effectively is lower by a factor of 12. Therefore in this case de- 
coupling is expected at a velocity va = 2 10 km/s and at a radius 
r — 1 . 1 3 R t . As the decoupling occurs near the stellar surface at a 
velocity that is lower by about a factor of five than the surface es- 
cape velocity v esc = 1 136 km/s, the passive plasma cannot be ex- 
pelled from the star. For metal abundances reduced by more than 
a factor of five with predicted mass-loss rates M $ lO~ 13 M /yr, 
the decoupling is even expected in the subsonic region. 

If the stellar parameters and the calcuated mass-loss rates 
are inserted into Eq. (50), with a solar number ratio of metals to 
hydrogen and helium y « 10~ 3 , it can be seen that this criterion 
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r eff [K] cr e [cm 2 g~'] Zi Z 2 
40000 0.33 1.08 2.95 

50000 0.34 1.09 3.51 

Table 2. Adopted electron scattering opacities and mean charges 
Zi (hydrogen and helium) and Z2 (metals). 



for ion decoupling leads to similar conclusions. The values of 
Vmax/Vesc are given in the last column of Table (1). It is v max > v ro 
only for solar metallicity. For the cases with reduced metallicity, 
it is v max < Vk,, although criterion (50) overestimates v max . Thus 
it follows that, at least for the cases with reduced metallicity, 
no coupled wind can exist. This could only be possible, if the 
terminal velocities were significantly lower than predicted from 
the present calculations. 

Krticka & Kubat (2000) suggest that in some cases decou- 
pling could be avoided if the wind switches to a solution with an 
abrupt lowering of the velocity gradient. However, the physical 
relevance of such solutions has been questioned by Owocki & 
Puis ( 120021 and Krticka & Kubat (120021) , who argue that they 
are unstable. The time-dependent hydrodynamical calculations 
of Votruba et al. ((2007 ) show that the ions decouple and are ac- 
celerated to high velocities, whereas the passive plasma deceler- 
ates. 

Feldmeier & Nikutta (2006) investigated the effect of ra- 
diative coupling between distant locations in winds with non- 
monotonic velocity laws and allowed for wind deceleration. 
They find that this could lead to a lowering of the terminal veloc- 
ity, which may be below the surface escape velocity by a factor 
of the order two to three. However, the existence of chemically 
homogeneous winds with M » lO~ 13 M /yr, which have been 
assumed e.g. by Unglaub & Bues (2001) to explain the helium 
abundances in sdB stars, would require the existence of winds 
with v'oo a 0.1v esc . Such a strong lowering of the terminal veloc- 
ity seems to be unlikely. As discussed in Sect. 6 from the various 
simplifications used in the present calculations, it should instead 
be expected that the terminal velocities are underestimated. Then 
the mass-loss rates required for the existence of a coupled wind 
would be higher. 

4. Results for r eff - 40000 and 50000 K 

This section presents the mass-loss predictions as a function of 
the surface gravity for T eff = 40000 K and 50000 K and for 
M» = O.5M (see Fig. 5). The radiative acceleration is ob- 
tained with the force multipliers from Kudritzki (2002) for sev- 
eral metallicities. The adopted values of n e /W are shown in the 
lower panels, which were obtained according to the iteration pro- 
cedure described in Sect. 2.4.1. For each wind model, it was 
checked whether decoupling of metals from hydrogen and he- 
lium is expected according to criterion (41); the adopted mean 
charges of these constituents are given in Table 2. The results of 
the calculations for which the actual one component description 
of the wind may be sufficient, are represented by solid lines. The 
surface gravities for which wind models were calculated range 
from values where the star would be near the Eddington limit 
(r e « 0.8) to values for which no wind solution exists. 

For r e ff = 40000 K, Z/Z = 1, coupled winds can exist up 
to \ogg = 5.8 with a mass-loss rate of M = 4.4 * lO~ 12 M /yr. 
For Z/Z = 0.2, this is only the case for logg < 5.3. Due to the 
higher luminosities for T e ff = 50000 K, the mass-loss rates for 
given log g and Z are higher than for = 40000 K. However, 
even for solar metallicity, chemically homogeneous winds can 
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Fig. 6. Predicted mass-loss rates (in M /yr) as a function of sur- 
face gravity for T eft - = 35000, 30000, 25000 K and Z/Z = 1, 
1/3, 0.1 with M, = O.5M . Dashed lines indicate decoupling of 
metals in the corresponding wind model. The upper and lower 
dotted lines in each figure represent the results of Vink & Cassisi 
(120021 ) for Z/Z = 1 and 0.1, respectively. 



exist for logg < 6.2 alone. Thus for DA white dwarfs in the 
considered range of effective temperatures, their existence can 
be excluded. For logg > 7.0 no wind solution exists at all, be- 
cause the radiative acceleration is too low even for wind opti- 
cal depth parameters t -> 0. For r eff = 50000 K, Z/Z = 0.2, 
and Z/Z = 0.01, we expect the existence of coupled winds for 
logg < 5.5 and for \ogg < 4.6, respectively. 

The results show an increasing dependence of the mass-loss 
rates on metallicity to higher gravities. For T e ff = 40000 K, 
\ogg = 4.0, and Z/Z = 1, a mass-loss rate of M = 1.3 * 
10~ 8 M o /yr is predicted. For Z/Z = 0.2 and the same stellar pa- 
rameters, it is M = 5.0* 10~ 9 M o /yr. From a comparison of these 
two values, a dependence of the mass-loss rate on the metallicity 
according to M ~ Z 58 may be suggested. For \ogg = 5.0 and 
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the same effective temperature the mass-loss rates for the two 
metallicities are 1.3 * 10~ 10 and 3.7 * lO~ n M /yr, respectively. 
From these values it would follow that M ~ Z 80 . This shows 
that, for weak winds, the Z-dependence cannot be represented 
by a power law of the form M ~ Z e with constant exponent e. 
For the present example, e would increase from e ss 0.6 for the 
lowest gravities up to values e > 1 for \ogg ^5.5. This effect 
can still be seen more clearly for the case T e g = 50000 K in Fig. 
5, for which results for Z/Z = 0.01 are also shown. This break- 
down of the usual scaling relations for sufficiently weak winds 
has already been found by Kudritzki (2002 ) in his calculations 
for stars with extremely low metallicities. 

From the middle panels of Fig. 5 it can be seen that the ter- 
minal velocities of the present wind models have a similar or- 
der of magnitude to the surface escape velocities. This is ex- 
pected from the original version of the CAK theory, as used in 
the present paper. According to the improved version of the the- 
ory they should be larger (see Sect. 6.1). With M(t) = kr a , 

it follows from the original theory that v M = J^z^ (see e.g. 

Lamers & Cassinelli 1 1999t . For a between about 0.5 and 0.7, 
ratios of 1.0 55 Vco/v^ & 1.5 should be expected. As for small 
wind optical depth parameters the log M(t) - log t relation fiat- 
tens and its slope -a finally approaches zero, the predicted ratio 
Voo/vesc ma y b e lower for weak winds, however. The sudden de- 
crease of Vco/v'esc and in the mass-loss rates, which occur in some 
cases near the highest gravities for which a wind solution still ex- 
ists (e.g. for T e ff = 40000K, Z/Z = 1 near logg = 6.3), are a 
consequence of the restriction of Kudritzki's force multipliers to 
a value M max as described in Sect. 2.4.1. The results for lower 
gravities, however, are not affected by this restriction. 

The results of this section show that the existence of coupled 
winds requires mass-loss rates at least of the order 10~ n M o /yr. 
For stars with M, = O.5M , effective temperatures in the range 
40000 K < T eff < 50000 K, and solar metallicity, this is possible 
only for surface gravities logg 6.0. Stars with stellar parame- 
ters within these ranges may be pre-white dwarfs, which are in 
an evolutionary stage between the asymptotic giant branch or the 
EHB and the white dwarf cooling sequence. Thus, if the metal- 
licity is not too far below the solar value, the mass-loss rates 
should be high enough to prevent the effect of diffusion. 

In no case can chemically homogeneous winds with M £ 
lO _12 M /yr exist, for which the effect of diffusion would be- 
come effective according to e.g. Unglaub & Bues (2001). That 
this is hardly possible can already be seen from Eq. (50). If 
e.g. for r eff = 40000K, logg = 6.0 and y = 10~ 3 (which ap- 
proximately corresponds to solar metallicty) a mass-loss rate 
M = 10~ 12 M G /yr is inserted; then with Z 2 » 3.0, it follows that 
Vmax/vesc = 1- Thus in winds with terminal velocities ^ v esc , 
which should be expected at least according to the improved ver- 
sion of the CAK theory (see Sect. 6), it is v max < v^,, so the 
constituents cannot be coupled throughout the wind. 

5. Results for sdB stars 

In Fig. 6 for T eB = 35000 K, 30000 K, 25000 K, Z/Z = 1, 
1 /3, 0.1, and for M» = O.5M , the predicted mass-loss rates are 
shown as a function of the surface gravity. These results have 
been obtained with force multipliers from our own calculations 
(see Sect. 2.4.2). The adopted values for the electron scattering 
opacity and for the mean charges of hydrogen and helium and 
for the metals are given in Table 3. 

For the typical range of surface gravities of sdB stars (5.0 ^ 
logg $ 6.0), the predicted mass loss rates are between about 
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Fig. 7. Force multipliers according to our calculations as a func- 
tion of the wind optical depth parameter for T e g = 30000 K and 
Z/Z = 1, 1/3, and 0.1. 
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Table 3. Adopted electron scattering opacities <x e , and mean 
charges Z\ of hydrogen and helium and Z2 of the metals, re- 
spectively. 



lO" lo M /yr and zero. The results strongly depend on surface 
gravity and metallicity. At least for such cases with M % 3 * 
lO~ 12 M /yr, decoupling of the metals from the passive plasma 
is expected. Thus the mass-loss rates may be high enough to 
prevent decoupling, but only in the most luminous sdB stars. For 
the more compact sdB stars from the one-component description 
of the wind mass-loss rates between about lO~ 12 M /yr and zero 
are predicted. However, the constituents cannot be coupled in 
such weak winds. 

Vink & Cassisi (2002) derived mass-loss rates from the re- 
quirement of a global momentum balance with the assumption 
of a yS-type velocity law with Voo = v esc - They calculated the 
radiative acceleration by means of a Monte Carlo simulation 
including about 10 5 of the strongest lines of the elements H- 
Zn with NLTE occupation numbers for the most important ele- 
ments. As can be seen from Fig. 6, especially for Z/Z = 1 and 
for 5.0 ^ logg S5 6.0, their results are in good agreement with 
the present ones. This may appear surprising, because Vink & 
Cassisis's calculations of the radiative force are clearly more so- 
phisticated. This point will be discussed in more detail in Sect. 
5.1 for the case r eff = 30000 K. 

The steep decrease in the mass-loss rates for Z/Z =0.1 and 
logg ^ 5.0 does not contradict Vink & Cassisi's predictions, 
as it may seem if their results for more luminous stars were ex- 
trapolated to higher gravities. Calculations of Vink (priv. comm.) 
for cases with logg = 5.5 and Z/Z = 0.1 have shown that no 
mass loss rate exists for which a global momentum balance can 
be fulfiled. Thus, if multicomponent effects are neglected, the 
mass-loss rate must approach zero for logg < 5.5. 



stars 





Fig. 8. Critical wind optical depth parameters as a function of the 
surface gravity for T eft - = 30000 K and Z/Z = 1, 1/3 and 0.1 . 

5.1. Detailed discussion forT eW = 30000 K 

In Fig. 7 for T eS = 30000 K and Z/Z = 1, 1/3, and 0.1, the 
force multipliers according to own calculations are plotted as a 
function of the wind optical depth parameter. The mass-loss rates 
shown in the middle panel of Fig. 6 were calculated with these 
force multipliers. The corresponding critical wind optical depth 
parameters (see Sect. 2.1) are shown in Fig. 8. It can be seen that, 
for the considered surface gravities 4.6 < logg < 6.2, they are 
in the range -5.5 > logf c ^ -8.0. As can be seen from Fig. 7, 
the slope -a of the log M (?) - log t relation is not constant in this 
range of wind optical depth parameters and a decreases to lower 
values of t. Thus, in a scaling law according to M ~ L~ (Puis et 
al. 120001 ), the exponent varies. An increasing dependence of the 
mass-loss rate on the luminosity to thinner winds is expected. 
In the present calculations for fixed T e g and M», this explains 
why the dependence of the mass-loss rate on the surface gravity 
increases to higher values of log g, as can be seen from the results 
in Fig. 6. 

In Fig. 9 for T eS = 30000 K and Z/Z = 1 and 0. 1, the contri- 
butions of the various elements taken into account in the calcula- 
tions presented in this section to the force multiplier are plotted 
as a function of the wind optical depth parameter. It can be seen 
that in all cases the major contribution is due to carbon for low 
values of f, which are relevant in the present calculations. From 
the present LTE assumption, it follows that 58 % of all carbon 
is C III and 42 % is C IV. The contribution of other ionization 
states of carbon is negligible. 

In Fig. 10 for Z/Z = 1, the contributions of the lines C III 
/1977.0A, /U175.7A and of C IV ,11549.1 A to the force multi- 
plier are shown as a function of t. At logf = -6.0, about 50 % 
of the total force multiplier is due to these three lines alone. For 
t — > this contribution increases to 77 %. The most important 
one of these three lines, especially for t — > 0, is CIII A971A. 

Now let us discuss the special example T e g = 30000 K, 
logg = 5.5. A coupled wind with a mass-loss rate M = 4.7 * 
lO~ 12 M /yr is predicted. The critical wind optical depth param- 
eter is log f c = -6.7. Here the contribution of these three carbon 
lines to the total force multiplier is 67 %. If all elements other 
than carbon were neglected, the derived mass-loss rate would be 
lower only by a factor of two (M = 2.5 * lO~ 12 M /yr). The effect 
of uncertainties in the CIII / CIV ionization equilibrium have a 
similar order of magnitude. If it were assumed that all carbon 
is in the ground state of C III, so that /1977A is the only line to 
significantly contribute to the radiative force, this again would 
lead to M — 2.5 * 10~ 12 M o /yr. If on the other hand it is assumed 
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Fig. 9. Force multipliers due to the elements H, He, C, N, and O 
for T e ff = 30000 K and for a) solar abundances, b) abundances 
of C, N, and O reduced by a factor of ten. 



that all carbon is in the ground state of C IV, this would reduce 
the mass loss rate by a factor of 3.6 to M — 1.3 * lO~ 12 M /yr. 
These results may explain the agreement between the present 
predictions and Vink & Cassisi's within a factor of about two. 
As the major contribution to the radiative force is due to a few 
lines of carbon, the number of other lines and elements taken 
into account is not essential. Secondly, the predicted mass-loss 
rates on the other hand always have the same order of magnitude, 
independent of the C III / CIV ionization equilibrium. If the ion- 
ization equilibrium changes, then the next ionization state takes 
over the contribution to the radiative force. 

If the metal abundances are reduced by a factor of ten, then 
the major contribution to the force multiplier is still due to car- 
bon (Fig. 9b). According to the present LTE calculations, the 
second largest contribution is due to hydrogen. For a wind op- 
tical depth parameter logf = -6.0, about 71 % of this contri- 
bution is due to the line La. Thus it essentially depends on the 
number density of the particles that are in the ground state of 
HI. The results shown in Fig. 9 were obtained with a number ra- 
tio of these particles to all hydrogen particles HI(n= 1)/(HI+HH) 
= 7.5 * 10~ 6 . If the HI/HII ionization equilibrium is shifted to 
HII, then, in contrast to the CIII/CIV ionization equilibrium, the 
next ionization state cannot take over. Thus a strong dependence 
on the ionization equilibrium is expected and the present results 
for hydrogen are questionable, because NLTE effects have not 
been taken into account. As the major contribution to the force 
multiplier for Z/Z = 0.1 is still due to carbon, however, the 
predicted mass-loss rates still agree approximately with the ones 
of Vink & Cassisi (2002). A possible overestimate of the rela- 
tive contribution of hydrogen to the force multiplier may be the 
reason that, for the cases with logg $5.0 shown in Fig. 6, the 
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Fig. 10. Force multipliers due to the lines CIII /1977.0 and 
,11 175 .7 A and CIV A1549. 1 A for T eff = 30000 K and solar abun- 
dances. The dotted line represents the total force multiplier. 



dependence of the mass-loss rates on the metallicity is somewhat 
less than in Vink & Cassisi's results. 

For r eff = 25000 K the situation is similar to r eff = 30000 K. 
The major contribution to the radiative force is due to carbon and 
the line C III ^977A. For T eS = 35000 K and wind optical depth 
parameters log t ^ -6.0 according to the present results, carbon 
is still the most important element. However, the contributions 
of N and O to the force multiplier are only slightly lower, by less 
than a factor of two. 



6. Discussion 

In Sects. 6.1 and 6.2 we discuss the effects of the omission of the 
finite disk correction factor, of neglecting changes in ionization 
in the wind, and of the shadowing of the flux by the photospheric 
lines. Possible consequences of the results for the surface com- 
position of the various types of chemically peculiar stars are then 
discussed in Sects. 6.3 and 6.4. 

6.1. Finite disk correction and changes in ionization 

The combined effect of the finite disk correction factor and 
changes in ionization may be estimated by comparing results 
from the present calculations with the theoretical predictions 
from other authors, who used the improved version of the CAK 
theory, which takes both effects into account. To eliminate the 
uncertainties of the force multipliers from our calculations, the 
values according to Kudritzki d2002l ) are used for this purpose. 
The inclusion of the finite disk correction alone should lead to 
lower mass-loss rates and higher terminal velocities by factors 
of 1.5 to four (Lamers & Cassinelli 1 1 9991 1 . 

Kudritzki (2002 ) calculated mass-loss rates for massive and 
luminous stars. Two examples of stellar parameters for the case 
Teff = 50000 K are given in Table 4. Ideally, our results for these 
stellar parameters should agree with Kudritzki's ones. As the 
finite disk correction and changes in ionization have been ne- 
glected in the present calculations, however, the results will usu- 
ally disagree. This can be seen from the comparison of both re- 
sults in Fig. 1 1 . In the last column of Table 4 the values of n e / W 
obtained iteratively as described in Sect. 2.4.1 are tabulated. In 
the present calculations these values are assumed to be constant 
throughout the wind. 

From the comparison it can be seen that our mass-loss rates 
in all cases are higher than Kudritzki's. In the various cases the 
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Table 4. Stellar parameters and adopted values of n e /W for 
the comparison with the mass-loss predictions from Kudritzki 
(12002b . with r eff = 50000 K in all cases. 
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Table 5. Comparison of the mass-loss predictions for Z/Z = 1 
with the results from Pauldrach et al. ((2004) in the upper part of 
the table and Pauldrach et al. ( 1988 ) in the lower. For each set of 
stellar parameters the results from Pauldrach et al. are given in 
the lower line. 



discrepancies vary by a factor of 1.5 up to a factor of seven. 
The terminal velocities are lower than Kudritzki's in most cases. 
Only for extremely low metallicities, the predicted terminal ve- 
locities approximately agree. 

These results correspond to the expectations. The effect of 
ionization may be more or less strong, depending on which el- 
ements and which ionization states primarily contribute to the 
radiative force. In some cases the combined effect of finite disk 
correction and ionization may lead to greater discrepancies than 
would be expected from the finite disk correction alone. In other 
cases, both simplifications may almost compensate each other, 
so that the results agree approximately with results from the im- 
proved version of the CAK theory. 

In Table 5 mass-loss rates and terminal velocities obtained 
with the present computational method and force multipliers ac- 
cording to Kudritzki ( 2002 ) are compared with predictions from 
Pauldrach et al. (120041 1 and Pauldrach et al. (119881 ) for central 
stars of planetary nebulae with a solar composition. The adopted 
values of n e /W are given in the last column. In one case only are 
the mass-loss rates from our slightly lower (by a factor of 1 .2) 
than the ones of Pauldrach et al. ( 1988 ). In all other cases, they 
are higher up to a factor of three. The terminal velocities in all 
cases are lower than the ones from Pauldrach et al. by factors be- 
tween about 1.5 and three. This again agrees with the expected 
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Fig. 11. Comparison of the predicted mass-loss rates (top) and 
terminal velocities (bottom) as a function of metallicity accord- 
ing to the present computational method (solid lines) with the 
results from Kudritzki (2002; dashed lines) for log L/L Q = 7.03 
and 6.42, respectively, with r ei f = 50000 K. In the present calcu- 
lations, Kudritzki's force multipliers have been used as described 
in Sect. 2.4.1. 

tendency of the predicted mass-loss rates of the present paper to 
be too large, whereas the terminal velocities are too low. 

6.2. The effect of line shadowing 

That the major contribution to the radiative force is due to a few 
strong lines for small wind optical depth parameters (see Sect. 
5.1) has important consequences. In the present calculations the 
emergent monochromatic flux F v at the frequency of a line, from 
which the force multiplier is obtained (see Eq. ( 1 9)), corresponds 
to the continuum flux. In Fig. 12, for the line CIII /1977.0A and 
for the case = 30000 K , logg = 5.5, and solar metal abun- 
dances, this flux is compared to the emergent flux from a LTE 
model atmosphere for sdB stars as used for spectral analyses 
(Heber, priv. comm.). It can be seen that near the line centre 
this flux is lower by about a factor of 100 than the flux that has 
been assumed in the mass-loss calculations. Thus in the inner 
parts of the wind where the velocity and thus the Doppler shift 
is low, F v has clearly been overestimated. For larger Doppler 
shifts, both fluxes approach each other and approximately agree 
for AA ^ 3. 5 A. These Doppler shifts correspond to velocities 
v ^ 1000 km/s, which may be compared with the surface escape 
velocity v esc = 955 km/s for this example with logg = 5.5. Thus 
the flux assumed in the calculations is approximately correct for 
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Fig. 12. Physical flux (in erg cirT V 1 ) near the line CIII /1977.0A 
for r e ff = 30000 K, logg = 5.5 and solar metal abundances as a 
function of AA = A - 977.0 A. The solid line represents the flux 
from a model atmosphere (Heber, priv. comm.), the dotted line 
represents the flux used in the present calculations. 



the outer parts of the wind where v <: v esc ; however, it is clearly 
overestimated in the inner parts. 

According to the results of Babel d 19961 1 this effect of "line 
shadowing" should lead to lower mass-loss rates and higher ter- 
minal velocities in comparison to the usual theory that assumes a 
velocity-independent flux. Then terminal velocities of the order 
Vk, « 5v esc and mass-loss rates that are lower by a factor of the 
order ten are not unusual. Thus we expect that the present cal- 
culations overestimate the mass-loss rates. The agreement with 
the results of Vink & Cassisi (|2002| l does not contradict this, 
because the terminal velocity is a free parameter in their cal- 
culations. If a terminal velocity = 5v esc instead of v = v esc 
were assumed, then it follows from the dependence of M on 
the assumed value of derived in Vink et al. (2000) that the 
mass-loss rates may indeed be lower by a factor around ten. The 
consequence of lower mass-loss rates would be that decoupling 
of the metals occurs at lower gravities than predicted according 
to the present results. 

Qualitatively, the inclusion of the effect of line shadowing 
should have similar consequences to including the finite disk 
correction (lower mass-loss rates, higher terminal velocities). 
Quantitatively, however, it seems to be much stronger for the 
weak winds discussed in the present paper. According to Owocki 
& ud-Doula (2004) the finite disk correction reduces the radia- 
tive force in the inner parts of the wind only by a factor of about 
two, whereas the example shown in Fig. 12 shows that line shad- 
owing may reduce the radiative force due to a strong line up to 
two orders of magnitude. 

The results of Babel d 19961 1 predict increasing terminal ve- 
locities for increasing surface gravities and otherwise fixed stel- 
lar parameters, whereas the present results shown in Sect. 4 pre- 
dict a decreasing tendency for Vco. The effect of line shadowing 
probably also changes the dependencies of the mass-loss rates 
and of the terminal velocities on the metallicity. 

6.3. Consequences for the chemical composition 

The results presented in Sects. 4 and 5 have shown that, in hot 
white dwarfs and in the majority of sdB stars (see Sect. 7), no 
chemically homogeneous winds can exist, at least if the metal- 
licity is not higher than solar. Thus, if any mass-loss exists, it 
should change the surface composition. The only exception may 
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be that the decoupling of the metals from the passive plasma 
occurs at a radius, at which the velocity of hydrogen and he- 
lium already exceeds the local escape velocity. Then hydrogen 
and helium can still be expelled from the star. For the case 
where decoupling occurs before the local escape velocity has 
been reached, Porter & Skouza ( 1999 ) propose a periodic time- 
dependent scenario. The passive plasma decelerates after decou- 
pling and eventually stalls. Thus a shell of gas is generated that 
reaccretes to the star. In the course of time, this scenario should 
lead to a depletion of metals near the stellar surface. 

Another possible scenario is a pure metallic wind with hy- 
drostatic hydrogen and helium. For main sequence A stars, this 
has been discussed by Babel ( 1995). If the stellar atmosphere is 
approximately in a stationary state, then the metals must have an 
outward movement. This leads to an outward frictional force on 
hydrogen and helium. If these elements are in hydrostatic equi- 
librium (vi = 0) and if we demand that their partial pressure 
decreases in an outward direction (dp\jdr < 0), then it follows 
from the momentum equation (24) for hydrogen and helium that 
the outward frictional acceleration must be less than the gravita- 
tional one: 



'coll 



(51) 



Otherwise the outward frictional force on the passive plasma due 
to the outflowing metals would be so strong that it cannot be in 
hydrostatic equilibrium. As the density in the stellar atmosphere 
is much higher than in the wind, it is a reasonable assumption 
that the drift velocities of the metals are small in comparison to 
the thermal velocity. Then the collisional acceleration on hydro- 
gen and helium can be calculated from the linear approximation. 
For x <sc 1 it is G (x) » \~h^x, so that the momentum exchanged 

per unit volume and unit time via Coulomb collisions may be 
written as 



|Ag| = n x n 2 l n \v 2 - V\ 



with 
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The quantity n\n 2 l\ 2 corresponds to the resistance coefficient K sl 
derived by Burgers ( [19691 see his Eq. 24.14). With Eqs. (29), 
(52), (53), with Vi = 0, and the equation of continuity for the 
mean metal, it follows for the collisional acceleration on hydro- 
gen and helium that 
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Then condition (51) can be written as 



M 2 < 4nGM* 



m\m 2 



(54) 



(55) 



If the mass-loss rate of the metals fulfils this condition, then the 
collisional acceleration on the passive plasma never can exceed 
gravity, independent of how high the velocity of the mean metal 
is (because the linear approximation leads to an upper limit of 
g^jj). With M, = O.5M , mean charges Z x = 1.0, Z 2 = 3.0, 
mean masses m\ = 1 3m„, m 2 - 14.6m p , and with In A = 12.0, 
it follows for a temperature T = 35000K that 

M 2 < 3 * 10- 16 ^ . 

yr 



Thus hydrogen and helium may be in hydrostatic equilibrium if 
the mass-loss rate of the metals is lower than about lO~ 16 M /yr. 
On the other hand, we know from the previous results that 
the existence of a coupled wind requires a total mass-loss rate 
greater than at least lO~ 12 M /yr. For an approximately solar 
mass fraction of the metals of the order of 10~ 2 , this means 
that a coupled wind requires a mass-loss rate of the metals alone 
M 2 > lO~ 14 M /yr. If, however, M 2 were somewhere in between 
these limits, e.g. M 2 w 10~ 15 M o /yr, then neither of the two cases 
seems to be possible. This value is too low for the existence of a 
coupled wind and too high for the case with hydrostatic passive 
plasma. Then hydrogen and helium may either fall back onto 
the star or be expelled, depending on at which radius decoupling 
occurs. Which of the various scenarios is the appropriate one 
essentially depends on the flow of the metals. In the following 
we discuss the case with hydrostatic hydrogen and helium and 
with a pure metallic wind. If the metals are trace elements and 
the densities in the wind are sufficiently low, then the mass-loss 
rates of the various metals should be independent of each other. 

If the mass-loss rate of a metal is greater than zero and if the 
atmosphere is approximately in a stationary state, the metal must 
have an outward velocity. If the effect of concentration gradients 
is weak, then the radiative force must not only balance gravi- 
tational settling, but also must compensate the inward frictional 
force on the metal due to collisions with the passive plasma. Due 
to saturation effects, the radiative force on a metal increases with 
decreasing abundance. Consequently, the abundances of these 
metals with non-zero mass-loss rates tend to be lower than the 
abundances obtained from diffusion calculations that assume an 
equilibrium between the downward gravitational force and the 
upward radiative force. Only for these metals, for which the 
mass-loss rate is sufficiently low, should the equilibrium abun- 
dances agree with the ones derived from spectral analyses. 

Seaton ( 1996 1999 ) presented the results of time-dependent 
diffusion calculations for iron group elements in the stellar en- 
velopes of HgMn stars and allowed for an outflow of these el- 
ements at the outer boundary. The abundances are predicted to 
vary with time. It may in principle happen that two stars with 
similar stellar parameters and surface compositions have differ- 
ent internal compositions. In the absence of strong concentration 
gradients, a maximum possible value for the outward flow of a 
metal can be derived for each depth, if the radiative acceleration 
is known as a function of the abundance. This maximum of the 
flow is small in regions where the metal has a noble gas con- 
figuration and thus the radiative acceleration is low. In Seaton 
(119991 ), these regions are referred to as barriers. How the surface 
abundance of a metal changes with time should depend essen- 
tially on its mass-loss rate and on the location of these barriers 
in the stellar envelope. 

The abundance patterns of the various types of chemically 
peculiar stars discussed in the present paper qualitatively agree 
with these expectations, although exceptions exist. For sdB stars, 
equilibrium abundances have been predicted and compared with 
measured ones e.g. by Bergeron et al. (119881 ), Charpinet et al. 
([T997j ), Ohl et al. (2000), and Behara & Jeffery (|2007l ). For 
some elements (e.g. Fe and N) there is generally good agree- 
ment, whereas for others (e.g. Si) the predicted values may be 
too high by several orders of magnitude. According to Chayer et 
al. (|2006 ), the measured abundances of Ge, Zr, and Pb are some- 
where between the solar value and the predicted equilibrium 
abundances, which for these elements exceed the solar value by 
about two orders of magnitude. The wide spread of the measured 
abundances may be an indication of a time-dependent scenario. 
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For hot hydrogen-rich white dwarfs, metal abundances have 
been measured or predicted e.g. by Barstow et al. ((2003 2005 ), 
Chayer et al. dl995al [I995b] l, Good et al. ([20 051), Schuh et al. 
( 120031 ), Dobbie et al. ( T2005i ), and Vennes et al. j20031l2006| ). For 
some elements (Fe and O), the agreement is close in many cases, 
for others (e.g. C, N, and Ni) the abundances predicted from 
equilibrium diffusion theory tend to be too large. An exception 
is silicon, for which the measured abundances are usually larger 
than the predicted ones. Thus in hot hydrogen-rich white dwarfs 
and in some helium-rich ones (Dreizler |1999l ) the situation is in 
some respect similar to sdB stars. The predicted abundances tend 
to be larger than the measured ones. 

Another example is found in the HgMn stars. The manganese 
abundances in the stars analysed by Jomaron et al. ( 1999) in all 
cases are lower than predicted from the equilibrium diffusion 
calculations of Alecian & Michaud d 19811 ) and greater than the 
solar values. Nitrogen is typically depleted by at least a factor of 
100, which is more than predicted from diffusion theory (Roby 
et al. I19991 ). The abundances of neon scatter from slight defi- 
ciencies to underabundances by one order of magnitude or even 
more (Dworetsky & Budaj 2000). This may point to some time- 
dependent process. The results of Seaton ( 1996 119991 ) for var- 
ious iron group elements show that the presence of an outflow 
at the outer boundary in general leads to better agreement. An 
exception to the expected tendency is mercury, for which abun- 
dances greater than predicted from equilibrium diffusion theory 
have been detected (Proffitt et al., 119991 ). However, at least in 
some HgMn stars, elements seem to be distributed inhomoge- 
neously over the surface (Hubrig et al. 2006 ), which complicates 
the problem. 

6.4. Selective winds or turbulence? 

As explained in the preceeding section, in the presence of pure 
metallic winds the abundances of the metals tend to be lower 
than predicted from equilibrium diffusion calculations. This 
should be so for overabundant, as well as for deficient, metals. 
An alternative explanation of the discrepancies between mea- 
sured and predicted abundances may be the presence of turbu- 
lence, which could stem from mixing processes like convection 
or stellar rotation. Then, however, the abundances of the defi- 
cient metals should be larger than predicted from equilibrium 
diffusion calculations (Vauclair et al. 1978), because turbulence 
tends to level out concentration gradients and reduces the effect 
of gravitational settling. Additional mixing due to turbulence 
outside of convection zones has been assumed by Richer et al. 
(2000) in their calculations for AmFm stars (which represent the 
continuation of the HgMn phenomenon at lower effective tem- 
peratures) to reduce the amplitude of the predicted abundance 
anomalies to a level that agrees with observational results. 

In hot hydrogen-rich white dwarfs and sdB stars, however, 
the abundances not only of enriched metals are lower than pre- 
dicted, but also those of deficient metals. In particular, the effect 
of turbulence alone can hardly explain why elements that should 
be overabundant according to diffusion theory are in fact defi- 
cient, as has been found for argon in two of the hottest known 
DA white dwarfs (Werner et al. 2007). Thus the scenario with 
selective winds seems to be the most likely one for these stars. 

The only element that is generally less deficient than pre- 
dicted from equilibrium diffusion calculations, is helium. This 
is so for these hot white dwarfs, in which helium is detectable 
(Vennes et al. [19881) sdB (Michaud et al. [19891 , and HgMn 
stars (Michaud et al. |19791 ). According to the results of Krticka 
(2006 ) for peculiar B stars, this can hardly be explained within 
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Fig. 13. Lines in the T e g - log g diagram above which chemi- 
cally homogeneous winds may exist for Z/Z = 0.1, 1/3 and 
1, respectively. Squares and circles represent the sdB stars anal- 
ysed by Maxted et al. (2001) and Lisker et al. (2005), respec- 
tively. Filled symbols represent the sdB's with peculiar H„ line 
profiles, which may indicate the presence of a weak wind. 

the framework of selective winds. Thus the helium deficiencies 
should be due to gravitational settling. It is possible that com- 
pletely undisturbed stellar atmospheres hardly exist and some 
turbulence is always present. This may explain that helium sinks 
more slowly than expected. However, it is still unclear what the 
origin of this turbulence is exactly. 

7. Summary and conclusions 

In the T e ff -log ^-diagram of Fig. 13 for 25000 K < T e ff < 
40000 K and for Z/Z = 1, 1/3 and 0.1, the lines are shown 
above which chemically homogeneous winds may exist accord- 
ing to the present results. The samples of sdB stars analysed 
by Maxted et al. d2001| ) and Lisker et al. ( 120051 ) shown in the 
figure have been analysed for peculiar H a line profiles, which 
may be interpreted as wind signatures according to Heber et al. 
(2003) and Vink (2004). As can be seen in the figure, such sig- 
natures have only been detected in the most luminous sdB stars. 
According to the calculations of stellar evolution by Dorman et 
al. ([1993 ), these ones are in a post-EHB stage of evolution. 

In the figure it can be seen that the majority of sdB stars 
populate a region in the T e g - log g diagram that is just below 
the line above which for Z/Z = 1 coupled winds may exist. If 
multicomponent effects are neglected, the results predict weak 
winds with mass-loss rates M ^ lO _12 M /yr. However, for such 
weak winds the momentum exchange between the metals, on the 
one hand, and hydrogen and helium, on the other, is not effective 
enough, because the densities in the wind are too low. Possible 
winds should be selective winds that lead to additional changes 
in the surface composition, which have not yet been taken into 
account in the diffusion calculations with and without mass-loss. 

According to the present results coupled winds can exist for 
the most luminous sdB stars, if the metallicity is not too low. 
This would explain why only in these stars wind signatures have 
been detected. The predicted mass-loss rates are of the order 
of 10~ n to lO~ 1() M /yr. This agrees with the predictions of the 
mass-loss recipe of Vink & Cassisi (2002), nevertheless the exis- 
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tence of these winds is still uncertain from the theoretical point 
of view. In Vink & Cassisi's calculations, the terminal veloc- 
ity has been a free parameter, and their results are for = v esc . 
The assumption of higher terminal velocities would lead to lower 
mass-loss rates. In the present calculations, the finite disk correc- 
tion and the shadowing of the flux by the photospheric lines have 
not been taken into account. Especially the latter effect seems to 
be very important. As discussed in Sect. 6, both simplifications 
should lead to the mass-loss rates being overestimated and the 
terminal velocities being underestimated. Then for a star with 
given stellar parameters and chemical composition, the densi- 
ties in the wind would be lower. Consequently decoupling of the 
metals from hydrogen and helium would occur at lower gravities 
than predicted, and the lines shown in Fig. 13 would be shifted 
to lower gravities. Frictional heating may have a similar effect. 
In addition, the presence of density inhomogenities in the wind, 
which are not taken into account in the present calculations, may 
possibly change the mass-loss rates. This effect of "clumping" 
is still investigated (Oskinova et al. 120071 1. For luminous stars 
there is some evidence that the mass-loss rates may need to be 
revised downwards (e.g. Bouret et al.|2005 1 Fullerton et al. 2006 1 
Martins et al. |20031 Puis et al. |2006t . 

To confirm the existence of coupled winds in the most lu- 
minous sdB stars, more sophisticated calculations are necessary. 
Because the sdB stars have a variety of abundances, it is impor- 
tant to show which elements contribute primarily to the radia- 
tive force. As explained by Puis et al. (2000), this should be the 
most abundant elements with the strongest lines in weak winds. 
According to the present results, the most important element is 
carbon for r e jf = 25000 and 30000 K. For higher effective tem- 
peratures the situation is less clear. 

If during the post-EHB evolution indeed winds with mass- 
loss rates M ^ lO~ n M /yr were initiated, the outer layers of 
the star will be removed rapidly. This should lead to the abun- 
dance anomalies gradually decreasing, so elements that are not 
detectable in sdB stars within the EHB band may reappear on 
the surface during the post-EHB evolution. This may explain the 
result of Edelmann et al. (2006) that, in those sdB stars with 
r eff £ 32000 K near the EHB band, no Si, Mg, and Al have been 
detected, whereas these elements are present in their more lumi- 
nous counterparts with similar effective temperatures, which are 
in a post-EHB stage of evolution. 

The results confirm the present picture of white-dwarf chem- 
ical evolution. The variety of compositions from hydrogen-rich 
to helium-carbon-oxygen-rich observed in pre-white dwarfs is 
due to the dredge up of processed material (Werner & Herwig 
2006). Mass-loss rates between about 10~ n and lO~ 6 M /yr are 
predicted for pre-white dwarfs with M, = O.5M , 30000 K < 
T eff < 50000 K, and 3.4 < logg <, 5.0, dependent on stel- 
lar parameters and metallicities. These mass-loss rates are large 
enough to prevent diffusion, at least if the metallicity is not re- 
duced by more than about a factor of ten in comparison to the 
solar value. No wind solution exists for white dwarfs on the cool- 
ing sequence with similar effective temperatures. Even for solar 
metallicity, the maximum possible radiative acceleration is too 
low by about one order of magnitude. Thus somewhere during 
the evolution with T e g > 50000 K, the mass-loss rates must de- 
crease to very low values so that chemically homogeneous winds 
can no longer exist. The possibility that metallic winds still exist 
may explain the result of e.g. Good et al. (2005 ) that neither dif- 
fusion calculations that assume chemically homogeneous winds 
nor diffusion calculations that assume the absence of any mass- 
loss can explain the measured abundances in DA and DAO white 
dwarfs very well. 



As discussed in Sect. 6.3, the existence of pure metallic 
winds with mass-loss rates M ^ lO~ 16 M /yr and with hydro- 
static hydrogen and helium may be a promising scenario to ex- 
plain the abundance anomalies of the metals in the various types 
of chemically peculiar stars considered in this paper. Then, how- 
ever, the problem remains, why the abundances of helium in gen- 
eral are larger than predicted from diffusion calculations, which 
assume an equilibrium between gravitational settling and radia- 
tive levitation. Possibly the stellar atmospheres are not com- 
pletely undisturbed after all. 

Hot (pre-) white dwarfs already have a variety of helium (and 
metal-) abundances before the onset of gravitational settling is 
expected. A dependence on the star's history may also play some 
role in sdB stars. According to the analyses of Edelmann et al. 
(2003 ) and Lisker et al. (2005) two distinct sequences of sdB 
stars seem to exist, which are characterised by an offset in the 
helium abundances. This phenomenon can hardly be explained 
with a single atmospheric effect. 
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